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Preface 

T became  interested  in  the  concent  of  the  aeromaneuvering 
orbital  transfer  vehicle  while  doing  research  on  Space  .Shuttle 
unner  stares.  The  potential  for  increased  eanability  over  purely 
propulsive  OTVs  is  so  much  greater  as  to  warrant  a detailed  inves- 
tigation. 

This  thesis  is  an  extension  of  work  performed  by  the  Tockheed 
Missiles  and  ^pace  Company  under  NAPA  contracts  NAP8-285R6  and 
NAP8-31^c?2t  and  is  based  on  the  recommendations  for  future  study 
contained  in  these  reports.  T felt  qualified  to  examine  the  problems 
of  navigational  accuracy  during  the  aeromaneuvering  phase  of  flight 
since  T had  completed  the  course  sequences  in  Navigation  and  in 
Guidance  and  Control  at  AFTT, 

T would  like  to  express  mv  appreciation  to  my  thesis  advisor, 

Ma.i  Richard  M.  Potter,  for  his  advice  and  critical  review  of  my 
work  during  the  preparation  of  this  report,  t would  also  like  to 
thank  Mr  J.P.  Hethcoat  of  the  Ms*. shall  Space  Flight  Center  for  his 
assistance  during  my  research  for  this  study. 


Robert  -I.  Chambers 
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Nomenclature 


The  MKS  system  is  used  in  this  study. 


AMOOS  - Aeromaneuverinft  Orbit- to- Orbit  Shuttle 
Ar  - aerodynamic  reference  area  * 15.69  m^ 

Cd  - coefficient  of  draft  ■ 2.5125 
- coefficient  of  lift  « 1.3457 


gain  coefficients 


D - draft 

GEO  - fteosynchronous  equatorial  orbit 


h - altitude 


h0  - reference  altitude  * 7.315  km 
L - lift 


LEO  - low  earth  orbit 
m - mass  of  vehicle  ■ 6804  kft 


r - inertial  position  vector 
u - unit  vector  or  control  variable 


V - inertial  velocity 


Vr  - relative  velocity 

Ve  - inertial  velocity  of  atmosphere 

R - bank  anftle 

Rc  - commanded  bank  anftle 

R0  - nominal  bank  anftle  - 90  deftrees 

y - fliftht  path  anftle 


usatsi£Jt>xi &jA  ■. 


• ^ i '•  **  - *>  * 


6 - see  Fig.  1 

u - gravitational  parameter 

p - atmospheric  density 

pr  - atmospheric  density  at  sea  level 

y - see  Fig.  1 
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Abstract 

A preliminary  investigation  was  conducted  to  determine  the 
navigational,  accuracy  required  by  the  Aeromaneuveriny  Orbit-to- 
Orbit  Shuttle  (AMOOS)  duriny  the  atmospheric  phase  of  fliyht,  The 
guidance  scheme,  which  is  the  same  as  the  one  developed  by  Lockheed 
in  the  original  AMOOS  study,  uses  the  parameters  of  velocity,  fliyht 
path  anrle,  and  density  altitude  to  correct  to  a nominal  trajectory. 
Density  altitude  is  obtained  from  atmospheric  density,  which  5s 
calculated  from  vehicle  acceleration.  Errors  in  velocity,  flight 
path  anrle,  and  acceleration  were  introduced  into  a three-deyree 
of  freedom  computer  simulation  of  the  vehicle  trajectory  usiny 
bank  anrle  commands  generated  by  the  yuidance  equations.  Deviations 
from  standard  atmospheric  density  were  taken  into  account. 

The  maximum  errors  allowable  that  still  permitted  the  vehicle 
to  attain  its  phasing  orbim  aooyee  altitude  of  720  km  (±100  km)  were 
determined.  Three  types  of  error  were  investigated : constant, 

sinusoidal,  and  random.  For  time- vary! ny  errors  the  frequency 
dependence  was  examined,  as  was  cross- couplinr  of  errors. 

It  was  concluded  that  the  ruidance  scheme  can  tolerate  fairly 
larye  errors  and  still  ruide  to  an  acceptable  apoyee  altitude.  The 
naviyational  accuracy  required  for  the  atmospheric  phase  of  flirht 
is  within  the  capability  of  present  day  astrionics. 
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NAVIGATIONAL  ACCURACY  REQUIREMENTS  OF 
AEROMANEUVERING  SPACE  VEHICLES 


I.  Int redact ion 


Background 

There  is  considerable  interest  in  geosynchronous  equatorial 
orbits  for  a variety  of  civil  and  military  applications,  such 
as  communication  links  and  surveillance.  Unfortunately,  it  will 
be  beyond  the  capability  of  the  Space  Shuttle  to  reach  these  orbits. 
The  Orbiter  vehicle  is  limited  to  an  altitude  of  300  km,  and  an 
upper  stage  is  required  to  reach  geosynchronous  altitude  of  35*800  km. 
An  upper  stage  vehicle  that  has  been  proposed  for  the  round  trip 
from  low  earth  orbit  (LEO)  to  geosynchronous  equatorial  orbit  (GEO) 
is  the  Space  'Pug.  However,  because  of  the  large  velocity  change  (AV) 
required  for  the  round-trip  mission,  about  8740  m/sec,  the  large 
propellant  mass  fraction  doer  not  allow  the  vehicle  to  carry  much 
payload. 

An  attractive  alternative  thst  has  been  presented  by  Lockheed 
is  the  Aeromaneuverinr  Orbit- to- Orbit  Shuttle  (AM00S)  (Ref  1,  Ref  2). 
The  advantage  of  this  concept  over  the  purely  propulsive  Space  Tug 
is  that  part  of  the  AV  required  for  slowing  and  returning  to  LEO  is 
provided  by  making  a braking  pass  through  the  atmosphere.  After 
exiting  the  atmosphere,  a short  rocket  burn  at  apogee  will  raise 
the  vehicle  perigee  out  of  the  atmosphere  and  establish  a phasing 
orbit  for  rendezvous  with  the  Shuttle  Orbiter.  The  total  AV  required 
in  this  case  is  6413  m/sec,  reducing  the  fuel  load  considerably 
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and  allowing  for  two  to  three  tines  the  payload.  The  basic  idea 
is  that  for  a snail  increase  in  structural  strength  and  the 
addition  of  a thernal  protection  system,  a larre  savings  in  the 
amount  of  fuel  required  can  be  realized. 

One  of  the  problems  associated  with  the  AMOOS  concept  is  that 
a very  accurate  atmospheric  flight  path  is  required.  Tf  the  vehicle 
is  one  kilometer  low,  it  will  reenter,  and  ix’  it  is  one  kilometer 
too  hitfh,  it  will  miss  the  apogee  of  its  phasing  orbit  by  thousands 
of  kilometers. 

The  vehicle  guidance  problem  is  complicated  by  the  fact  that 
the  density  of  the  upper  atmosphere  varies  widely  at  hirh  altitudes. 
The  density  will  vary  from  +4C K to  -35^  from  its  nominal  value 
(Ref  1:67,  Ref  3^20).  The  AMOOS  vehicle  is  able  to  compensate  for 
off-nominal  conditions  by  modulating  its  lift  vector,  and  this 
aeromaneuverinr  ability  is  what  differentiates  this  proposal  from 
earlier  ballistic  trajectory  concepts. 

Purpose  of  the  Study 

One  of  the  areas  of  the  AMOOS  proposal  that  has  not  been 
investigated  so  far  is  the  accuracy  required  of  the  vehicle  nav- 
igation system  to  successfully  implement  the  ruidance  equations, 
and  this  is  the  subject  of  my  thesis.  T nave  taken  the  ruidance 
scheme  used  in  the  oxurinal  study  and  incorporated  it  into  a 
trajectory  simulation  in  which  navi ration  system  errors  can  be 
introduced. 

The  real  of  the  ruidance  system  is  to  place  the  vehicle  in  an 
orbit  with  an  aporee  of  720  km  it  100  km).  The  ability  of  the 
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guidance  system  to  accomplish  this  task  is  related  to  the  navigational 
errors.  The  navigation  parameters  in  the  guidance  calculations 
that  are  subject  to  errors  are  the  vehicle  velocity,  flight  path  angle, 
and  acceleration.  In  the  simulation,  errors  in  these  variables  are 
introduced  in  three  ways?  as  a constant  error  in  navigational 
information,  as  a sinusoidal  error,  and  as  a random  error.  The 
results  of  the  simulation  can  be  used  to  design  a navigation  system 
capable  of  supplying  velocity,  flight  path  angle,  and  acceleration 
information  within  the  required  accuracy. 

Previous  Investigations 

A previous  navigational  accuracy  study  has  been  accomplished 
for  the  aerobraking  return  from  0150  to  LEO  (Ref  4:33-45).  However, 
this  study  examined  a vehicle  that  followed  an  unguided  trajectory 
through  the  atmosphere  (it.  used  aerobraking,  but  not  aoromaneuvering). 
The  results  of  this  study  indicated  that  a one-pass  unguided  braking 
maneuver  was  not  feasible  using  forecast  shuttle-era  technology, 
due  to  the  high  accuracy  requirements  and  unpredictable  variations 
in  atmospheric  density. 

A preliminary  investigation  by  Lockheed  of  the  effects  of 
exoatmc.spheri c navigation  errors  has  determined  that  it  will  be 
necessary  to  position  the  AMOOP  vehicle  within  a corridor  of  + 3.5  km 
in  depth  nrior  to  entry  (Ref  2:73).  The  investigation  did  not  include 
the  effects  of  atmospheric  navigation  errors.  These  errors  are  the 
subject  of  this  study,  which  solely  concerns  itself  with  navigational 
errors  generated  within  the  atmosphere  and  their  effect  on  the  guidance 
equations.  Of  course,  the  two  studies  are  not  independent,  since  the 


same  astrionies  will  most  likely  be  used  in  both  cases.  It  is 
expected  that  the  astrionic  equipment  will  include  a combination 
of  the  following  sensors:  Inertial  measurement  unit  (IMIJ),  star 

tracker,  horizon  sensor,  landmark  tracker , and  laser  or  conventional 
radar. 

It  is  not  the  purpose  of  this  study  to  actually  design  a 
na vibration  system  and  select  the  equioment  used  (each  of  which 
would  have  its  own  characteristic  errors),  but  the  results  indicate 
to  what  tolerances  the  navigation  system  must  supply  information 
to  implement  the  guidance  scheme. 

Chapter  II  will  present  the  guidance  law  used  in  the  study. 
Chapter  III  will  be  a description  of  the  simulation,  and  Chapter  TV 
covers  the  introduction  of  navigational  errors.  Results  are  presented 
in  Chapter  V,  and  conclusions  are  in  Chapter  VI. 


Explanation  of  Method 

The  AMPOS  vehicle  is  constrained  to  fit  within  the  Shuttle 
Orbiter  payload  bay,  and  is  basically  a cylinder,  with  one  end 
truncated  into  a lifting  body  shape,  A detailed  description  of 
the  vehicle  is  included  in  Appendix  A.  For  the  purpose  of  deriving 
the  guidance  law,  it  should  be  noted  that  the  vehicle  has  a L/D 


ratio  of  from  .5  to  1.0  at  hypersonic  velocities. 

The  guidance  problem  can  be  simply  stated  a maneuvering 
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the  vehicle  so  as  to  leave  the  atmosnhero  with  .fast  enough  velocity 
to  hit  the  proper  apogee  of  the  phasing  mbit  (720  km).  This  is 
a different  problem  than  conventional  reentry  guidance,  which  is 
concerned  with  reducing  velocity  below  some  terminal  value,  with 
heating  and  dynamic  pressures  being  the  primary  constraints. 

The  guidance  method  that  has  been  proposed  for  AM00S  is 
controlling  to  a nominal  trajectory,  using  a linear  regulator. 

The  lift  vector  is  used  to  steer  to  the  nominal  trajectory  by  bank 
angle  modulation.  Rank  ancle,  the  orientation  of  lift,  is  used 
instead  of  ancle  of  attack  for  lift  control  because  the  vehicle  shape 
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its  inclination  from  zero  to  28.5  decrees  (the  orbital  inclination 
of  the  Shuttle  Orbiter,  and  the  latitude  of  Kennedy  Space  Center). 

A nominal  bank  anple  of  90  decrees  produces  a plane  change  capability 
of  approximately  7 decrees.  This  reduces  the  propulsive  requirements 
at  geosynchronous  altitude  for  plane  change  aV.  Changes  in  the 
bank  an^le  from  90  decrees  for  steering  corrections  will  reduce  this 
plane  change  capability  slightly. 

Fquations  of  Motion 

The  parameters  used  to  generate  bank  an?le  commands  are  inertial 
velocity,  flight  path  anple,  and  density,  which  is  a function  of 
altitude.  The  original  derivation  of  the  ruidanee  law,  from 
Reference  1 , is  presented  here  because  it  provides  the  basis  for 
the  remainder  of  the  study.  The  equations  of  motion,  in  the 
osculating  plane  of  the  orbit,  are  derived  as  follows  (Ref  lj,D-t  - 
D-21  ) : 


dV  dV  de 

Al  * 3„dt  + «xV  dt 

where  u„»  unit  vector  in  the  direction  of  V 
ux«  unit  vector  perpendicular  to  V 
with  the  inertial  reference  frame  as  shown  in  Fir.t. 
Summinr  the  forces  in  the  direction  of  V reives 


-D  - pm  sin  y 


m V 


(2) 


or  v - - - - -“a sin  y (3) 

m * 

Summing  the  components  perpendicular  to  V chives 
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L - ^am  cos  y ■ m V 9 
r* 


As  can  be  seen  from  Fig.  1 


■y  + y 


_\  _ V COS  Y 

v ■ — ? 


Combining  (5)  and  (6)  gives 


* . V 

0 - v — cos  v 

r 


Substituting  (7)  into  (4)  and  solving  for  y yields 


. i r v j*i 

V [ r*  r J 


cos  Y 


The  derivative  of  radial  distance  r is  given  by 


r » V sin  y * h 


If  the  states  of  the  system  are  taken  to  be  V,y,  and  h,  then 
equations  (3),  (8),  and  (9)  provide  the  state  equations. 

Although  the  lift  and  drag  force  directions  are  referenced  to 
inertial  velocity,  they  may  be  approximated  by  the  force  directions 
with  respect  to  the  vehicle  velocity  relative  to  the  atmosphere 
since  the  directions  of  relative  velocity  and  inertial  velocity  are 
almost  parallel  during  the  period  of  atmospheric  flight. 

Thus 


D a:  \ p V*  A_  C. 


(10) 


L ^ i p vi  Ar  C}  cos  R 


Vr  * V - vft 


(12) 


where  bank  ansle,  R,  is  defined  in  Fi*.  2. 

Control  of  the  vehicle  is  accomplished  by  rotating  the  lift 
vector  through  a bank  an^le  P,  which  changes  the  component  of  lift 
in  the  Yr  direction,  as  show  in  Fif?.  2.  The  relation  between  R 
and  the  vehicle  orientation  is  shown  in  Fig,  14,  Appendix  A.  The 
control  variable,  u,  is  defined  as  the  deviation  of  the  commanded 


value  of  C^y  from  some  nominal  value  of  bank,  i.e,, 

i,  - Cxy  - Ciy0  - Cx  ( cos  % - cos  fl0  ) 

where  RQ  «•  90  decrees  for  maximum  plane  change  caoability. 

Substituting  (10),  (11),  (12),  and  (13)  into  equations  (3) 
and  (8)  and  assuming  R - Rp  results  in  the  revised  equations  of 
motion? 


(n) 


SS  P (V  - Vefcd  Ap  -*S*stnv 

2mV  P (v  * ve^  Ar  fcl  sos  po  + - jpr 

f * V sin  v 


COS  Y 


(14) 


(15) 


(16) 


where  densitv  is  an  exponential  approximation  according  to 


P "Po  e 


-(h/h0) 


(17) 


These  equations  (with  u*0)  are  integrated  from  initial 
conditions  selected  to  result  in  an  anoree  altitude  of  720  km, 
Since  entry  altitude  i-  fixed  (at  120  km),  all  possible  initial 
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U - Cl  (V  - Vnotn)  + C2  (v  - Ynom)  + C3  (h  - hno!fl)  (IB) 

The  time- varying  coefficients  Cg,  and  C3  are  the  solution 
to  the  linear  reeqilator  problem  presented  in  the  next  section.  The 
values  of  Cp,  C^,  VnoyTlt  Ynom,  and  hn0Tn  can  be  precomputed  for 
discrete  time  intervals  and  stored  in  the  on-board  computer. 

Linear  Regulator  Problem 

The  developement  of  the  linear  remilator  problem  involves  a 
linearization  of  the  equations  of  motion  to  the  form  of 
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and 


(21) 


The  matrices  A and  B are  evaluated  along  the  nominal  trajectory, 
The  partial  derivatives  are  calculated  from  equations  (14),  (15),  and 
(16)  as  follows; 


|V  . . C v.  Cd  A, 


(22) 


av  m 


av  _ u 


av 


COS  Y 


(23) 


av  i 


f?  - - in  v?  cd  *r  §r  + r“  81,1  Y 


(24) 


dY 

av 


i 

mV 


[ 


V*  - 


-Vi 

2V 


p Ar  f Cl  oos  % + U | + j^^r-a  + ^ 


cos  Y (25) 


av 

*Y 


V - H - 

V V* 

m 


V 

r 


sin  Y 


(26) 


av 

ar 


2mV 


Vr  Ar 


cos  ft0  + u 


31. 

3r 


2U  + 
V r* 


?■] 


COS  Y 


(27) 


af 

av 


sin  Y 


(28) 


af 

aY 


V cos  Y 


(29) 
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f£  - 0 (33) 

ou 

and 

|P-  - - — (34) 

dr  h0 

Ilsinr  standard  optimal  control  techniques  (Ref  5:209),  a 
performance  index  is  chosen  to  minimize  the  erroi  in  phasing  orbit 
apogee.  This  index  of  performance  (which  is  to  be  minimized) 
consists  of  the  terminal  miss  which  relates  directly  to  the  square 
of  the  aoopee  error,  Plus  the  integral  of  the  square  of  the  control 
variable  u,  jvsr  the  nominal  flight  time.  This  last  term  tends  to 
limit  the  control,  hence  limiting  the  amount  of  propellant  required 
for  the  reaction  control  system. 
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The  solution  to  the  optimization  nroblem  is  given  by 


u - - R*1  BT  K 6x  - CT  6x 


whore  K Is  Pi von  by  the  Rioeati  differential  equation 


(3?) 


K - -At  K - K A + K R R-1  BT  K 


(33) 


which  is  integrated  backwards  through  time  from  its  final 
conditions,  K(tf)  - H. 

Weighting  Matrices 

The  matrix  H is  chosen  to  minimize  the  effect  of  the  final 
state,  5x(tf),  at  atmospheric  exit,  on  apogee  altitude  ha.  For 
small  variations,  &ha  can  be  related  to  the  final  state  by 


ih*  jv  + ih»  6v  + di»  8r 


t*tf 


(39) 


The  matrix  H is  determined  by 


6ha  - 6x^(tf)  H 6x(tf) 


(40) 


which  yields  the  elements  h^  of  H; 


h22”  OUl* 


h - *ha 
33  17 


u u ^ha  ^ha 

h10  - h_,  - — S — » 

12  21  av  a v 


(41  ) 


(42) 


(43) 


(44) 


1 , dha  dhn 

^ 3i  dv  dr 


and 

, « dha  Sha 

h23m  h32*  dy"  W 


(45) 

(46) 


Evaluating  the  partial  derivatives  at  the  nominal  conditions  at 
the  final  time  (V-7980  m/sec, Y *1.75  oeg,  and  r»6.49R,153  m)  results 
i.n  numerical  values  of 


dY 


74.73 


km 

deg 


?!% 

dr 


3.116  ® 

m 


The  function  R(t)  was  selected  to  prevent  the  peaking  of  the 
gains  at  perigee,  the  period  of  maximum  control  effectiveness. 

The  function  selected  is  shown  in  Fig,  3. 

The  feedback  gain  vector  derived  from  eq.  (37)  is  shown  in 
Fig.  4.  In  Figures  3 and  4,  the  point  at  which  time  t*»0  is  when 
the  vehicle  passes  through  95.4  km  upon  entry.  Above  this  altitude 
there  is  not  sufficient  atmosphere  to  permit  control  effectiveness. 

Tn  the  region  of  atmospheric  guidance,  the  vehicle  acceleration 
is  primarily  affected  by  the  drag  terms.  Therefore,  neglecting 
the  gravitational  terms,  and  for  small  flight  path  angles,  from 
eq.  (14) 
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Feedback  Gain  Vector  C for  AM 0G3  Guidance  (Ref  2:  D-12) 
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This  density  can  be  related  to  an  altitude  through  the  use  of  eq,(17) 


h - - h0  ln(  p/p0  ) 


(49) 


Thus,  from  measuring  vehicle  acceleration,  a density  altitude  can 
be  calculated.  This  is  compared  to  the  nominal  density  altitude  to 
determine  the  vehicle  altitude  error. 
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ITT.  Description  of  Simulation 


Purpose 

The  valid! tv  of  the  guidance  scheme  is  tested  by  a computer 
program  which  simulates  the  actual  vehicle  dynamics  with  guidance 
commands.  Perturbations  from  the  nominal  trajectory  are  introduced 
to  determine  how  well  the  guidance  law  will  oorreot  back  to  the 
nominal.  In  the  case  under  study,  perturbations  are  caused  by 
deviations  in  atmospheric  density  and  by  errors  in  the  calculated 
position  and  velocity  information  supplied  by  the  vehicle's 
navigation  system.  Atmospheric  density  fluctuations  are  known 
to  vary  from  -35$  to  4-0$  from  normal  in  the  region  flown  by  AMOOS 
(from  120  to  70  km),  and  are  easily  input  to  the  simulation. 
Navigational  errors  can  also  be  easily  input  to  the  simulation 
and  these  error  sources  are  described  in  the  next  chapter. 

Assumptions 

The  following  assumptions  have  been  made  rn  constructing  the 
vehicle  simulation: 

1.  The  assumptions  of  a spherical  rotating  earth  and  atmosphere. 

2.  The  atmospheric  density  can  be  modeled  as  an  exponential 
function  of  altitude. 

3.  The  vehicle  is  returning  from  GRO  via  a Hohmann  transfer  orbit 
with  a 21,5  decree  plane  change,  with  a targeted  vacuum  perigee  of 
71,1  km.  This  value  for  periree  is  the  result  of  adjusting  vacuum 
peripee  to  determine  the  initial  entry  conditions  that  result  in 
an  apopee  of  720  km. 
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k.  Vehicle  flies  constant  angle  of  attack  with  lift,  longitudinal 
vehicle  axis,  and  relative  wind  always  in  the  sime  plane.  The  orientation 
of  this  plane  (R)  is  the  only  vehicle  control  that  affects  translational 
dynamics , 

5.  The  actual  bank  angle  is  equal  to  the  commanded  bank  angle. 

It  was  originally  intended  to  simulate  an  autopilot  in  the  Guidance 
loop  to  implement  the  bank  commands,  and  to  account  for  the  vehicle* s 
inability  to  instantly  assume  the  commanded  value.  However,  from 
the  results  of  the  simulation,  the  guidance  commands  are  fairly 
continuous,  and  the  value  of  (P  - PQ)  is  small  (about  2 - 3 decrees). 

This  permits  the  vehicle  to  asstime  the  value  of  Pc  within  the  guidance 
cycle  time  (4  sec),  assuming  an  angular  acceleration  of  1 deg/sec  . 


Equations  of  Motion 

Since  it  was  desired  to  examine  the  three-dimensional  vehicle 
motion  rather  than  just  the  motion  in  the  osculating  orbital  plane, 
the  equations  of  motion  used  in  generating  the  guidance  law  could  not 
be  used.  A different  (although  consistent  1 set  of  equations 
was  used,  based  on  a inertial  geocentric-equatorial  (T.TK)  coordinate 
system.  The  T.TK  reference  frame  was  oriented  with  the  T direction  in 
the  direction  of  the  vacuum  neriree,  to  simplify  the  transformation 
from  neri focal  coordinates  to  T.TK  coordinates.  These  equations  are 
the  same  as  used  in  the  original  study  (Ref  lrA-3). 

The  vector  equation  of  motion  is: 
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where  I.  and  D are  the  unit  vectors  in  the  directions  of  lift  and  drag, 
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The  aerodynamic  coefficients  C]  and  and  the  vehicle 
specifications  of  mass,  m,  and  reference  area,  Ar,  are  taken  from 
the  design  of  Reference  2.  Relative  velocity,  Vr,  is  defined  by 
eq.  (12),  Bank  angle,  R,  comes  from  the  guidance  law  (eq.  (18)), 
and  the  density  altitude,  h,  is  computed  from  eqs.  (48)  and  (49). 
The  program  user  a Runge-Kutta  integration  algorithm  to  integrate 
th3  equations  of  motion  from  atmospheric  entry  to  exit  (120  km). 


Computer  Program 

A flow  chart  of  the  ma.ior  elements  of  the  computer  program  is 
shown  in  Fig.  5.  A listing  of  the  program  is  included  in  Appendix  B, 
The  program  first  generates  the  nominal  tra.iectory,  using  the  equations 
of  Chanter  TT.  The  vector  of  time- varyinm  guidance  coefficients  is 
input  as  data  points,  and  an  interpolating  table  look-up  routine  is 
used  to  reconstruct  the  curves  of  Fig.  4 and  select  values  for  a 
given  flight  time.  The  elements  V,  v,  and  r of  the  the  osculating 
orbit  at  r « 120  km  are  converted  to  inertial  position  and  velocity 
and  are  used  as  initial  conditions  for  the  integration  of  the  equations 
of  motion  given  in  the  previous  section. 

Simulated  navigation  errors  in  V,Y,  and  acceleration  are  input 
to  the  guidance  section,  which  then  generates  a bank  angle.  The 
integration  proceeds  until  the  vehicle  leaves  the  atmosphere 
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(when  It.  amain  reaches  120  km).  At  this  noint  the  apogee  of  the 
phasing  orbit  is  calculated,  assuming  two  body  dynamics. 

The  integration  time  step  is  1 second,  which  was  determined 
to  give  good  accuracv  without  resulting  in  excessive  computational 
time.  Mew  guidance  commands  are  generated  every  4 seconds,  and 
bank  angle  commands  are  limited  to  ± 10  degrees  about  the  nominal 
90  degree  bank  angle.  The  guidance  cycle  time  and  maximum  bank  angle 
represent  the  values  that  gave  the  best  results  in  the  original 
simulation  work  done  by  Lockheed.  A sample  trajectory  is  given  in 
Appendix  C. 
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T v • Navigation  System  Errors 

It  was  attempted  to  model  errors  in  a manner  that  would  be 
typical  of  the  outputs  from  the  navigation  system.  Accordingly, 
the  errors  were  input  to  the  simulation  as  a constant  (bias')  error, 
a sinusoidal  error,  and  as  a random  error.  Each  type  of  error  was  run 
with  varying  magnitudes  to  determine  its  relation  with  the  output, 
apogee  error.  Tn  addition,  each  error  in  velocity,  flight  path  angle, 
and  acceleration  was  considered  for  three  different  atmospheric 
density  profiles.  The  first  profile  was  the  nominal  density,  and 
the  second  and  third  profiles  were  the  worst-case  situations  of 
density  being  40^  greater  than  normal,  and  bein^  35^  less  than  normal. 
This  demonstrates  the  capability  of  the  guidance  system  to  correct 
for  off-nominal  atmospheric  conditions,  in  addition  to  navigation 
errors. 

The  sinusoidal  errors  were  run  at  varying  frequencies  to 
determine  if  there  was  any  frequency- sensitivity  in  the  guidance 
system.  The  range  of  frequencies  to  be  examined  was  determined  by 
the  guidance  cycle  time  and  total  time  of  flight.  The  high  frequency 
limit  was  chosen  to  give  a period  equal  to  the  guidance  cycle  time 
of  4 seconds.  Any  frequency  greater  than  this  would  appear  as  a 
random  error  to  the  simulation.  The  low  frequency  limit  was  selected 
to  give  a period  approximately  equal  to  the  time  of  atmospheric 
flight  (400  seconds).  A frequency  much  lower  than  this  would  appear 
to  be  almost  a constant  error. 

The  effects  of  random  errors  in  velocity  on  apogee  altitude  were 
examined  to  determine  if  there  was  any  significant  difference  between 


the  results  for  random  errors  and  sinusoidal  or  constant  errors. 


The  random  errors  were  produced  using  the  random-number  utility 
subprogram  of  the  computer  system.  A uniformly  distributed  random 
sequence  of  numbers  between  -1  and  +1  was  assigned  to  the  elements 
of  an  array  at  a discrete  interval.  This  interval,  which  could  be 
varied,  was  called  the  correlation  time.  The  elements  of  the  array 
between  these  intervals  were  filled  in  using  a linear  function 
between  the  random  end  points  of  the  interval.  The  elements  of  the 
array  were  multiplied  by  the  appropriate  error  magnitude  and  then 
each  element  corresponded  to  an  error  value  at  each  cycle  of  the 
guidance  equations.  The  use  of  the  correlation  time  provided  a means 
of  varying  the  rate  at  which  the  random  errors  were  input  to  the 
simulation. 

The  correlation  time  was  analogous  to  the  period  of  the 
sinusoidal  errors,  and  so  the  range  of  correlation  times  was  chosen 
from  4 to  400  seconds,  using  the  same  reasoning  as  was  used  to  determine 
the  frequency  of  '•inusoidal  errors. 

Errors  in  each  parameter  were  first  examined  independently, 
but  the  effects  of  cross-coupling  rf  errors  v*»re  also  examined. 

This  was  to  determine  if  the  principle  of  superposition  would  be 
applicable  to  the  error  sources. 

The  results  of  the  study  are  presented  in  the  next  chapter. 
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V.  Results 


Constant  Error 

The  nominal  tra.iectory  produced  by  eqs.  (14),  (IS),  and  (16) 
resulted  in  an  apogee  altitude  of  exactly  720  km.  However,  when  this 
tra.iectory  was  "flown”  in  the  three-dimensional  simulation  using  eqs. 
(SO),  (SI),  and  (52),  with  no  errors  input,  it  consistently  yielded 
an  apogee  altitude  of  V)0  km.  This  10  km  bias  probably  results  from 
the  different  methods  of  calculating  Vg  in  each  case.  In  the  case  of 
the  nominal  tra.iectory,  it  is  approximated,  since  the  orientation  of 
the  osculating  orbit  is  not  known,  while  in  the  3-D  simulation  it  is 
calculated  exactly.  This  bias  was  taken  into  account  when  determining 
the  results  of  the  navigational  errors. 

Errors  in  apogee  due  to  errors  in  measured  inertial  velocity 
exhibit  a linear  relationship,  as  shown  in  Fig.  6,  The  horizontal 
lines  show  the  permissible  variation  (+  100  km)  in  apogee  altitude, 
and  the  three  diagonal  lines  show  the  three  cases  of  atmospheric 
density  being  equal  to  nominal,  and  +40^  and  -35^  from  nominal. 

From  the  figure,  taking  into  account  the  10  km  bias  and  the 
possible  variation  in  atmospheric  densitv,  It  can  be  determined  that 
the  allowable  error  in  measuring  velocity  is  + 22  m/sec.  A .sensitivity 
coefficient  of  -4.25  km/m/sec  can  be  calculated  for  constant  errors 
(meani ng  that  a +1  m/sec  error  in  measuring  velocity  will  result  in 
an  apogee  error  of  -4,25  km). 

Figure  7 shows  the  effects  of  a constant  error  in  measured 
flight  path  angle.  This  time,  however,  the  effects  are  seen  to  be 
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Constant  Velocity  Errors 


non-linear.  The  maximum  allowable  errors,  again  allowing  for  the 
bias  error,  can  be  seen  to  be  -.0045  and  .0035  (-.26  deg  and  .20  deg). 

A linear  sensitivity  coefficient  cannot  be  determined  for  these  curves. 

Errors  in  measuring  acceleration  result  in  the  apogee  error 
shown  in  Fig.  8.  The  curves  indicate  that  acceleration  errors  must 
be  kept  between  ± .6  m/sec2.  The  curves  can  reasonably  be  approx- 
imated by  a straight  line  in  a small  region,  permitting  the  cal- 
culation of  an  approximate  sensitivity  coefficient.  For  the  region 
of  m/sec2  error,  the  sensitivity  coefficient  is  120  km/m/seo2. 

Sinusoidal  Error 

Sinusoidally- varying  errors  in  velocity,  flight  path  angle, 
and  acceleration  were  input  to  the  simulation  at  several  different 
amplitudes  for  each  of  the  three  atmospheric  density  profiles. 

Typical  results  are  presented  in  Figures  9,  10,  and  11.  These  are 
for  the  nominal  density  profile,  with  amplitude  of  the  errors 
being  respectively  20  m/sec,  .004  radians,  and  .4  m/sec2. 

The  figures  will  be  discussed  together,  since  they  all 
indicate  the  same  trend.  The  data  is  plotted  as  discrete  points 
because  it  did  not  lend  itself  to  fitting  a smooth  curve  between 
the  points.  Tt,  can  be  seen  that  at  frequencies  greater  than  .1  rad/sec 
the  apogee  deviations  are  small  compered  to  the  deviations  produced 
with  constant  errors  of  the  same  magnitude.  Below  .1  rad/sec 
the  deviations  increase  erratically  until  they  are  approximately 
the  same  as  for  constant  error.  This  is  most  likely  due  to  the 
guidance  scheme  being  more  sensitive  to  errors  at  one  point  in  the 
flight  oath  (probably  oeriree,  since  this  In  the  region  of  greatest 
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Sinusoidal  Flight  Path  Angle  Error  (amplitude  « .004  radians) 


'PfWlt’ifT'irfr-J'TTl 


wpwnr 


control  effectiveness).  The  high  frequency  errors  tend  to  cancel 
each  other  out,  but  at  the  lower  frequencies  there  is  a larger 
error  in  the  sensitive  region  that  is  not  cancelled  out.  The  net 
effect  is  that  the  system,  like  most  physical  systems,  acts  like  a 
low-pass  filter.  The  guidance  scheme  did  not  show  any  instability 
in  the  frequency  range  examined. 

The  results  using  the  off-nominal  density  conditions  (-3*5^, 
+40^)  were  the  same,  with  the  data  points  being  shifted  up  or  down 
as  in  the  constant  error  oases.  Varying  the  amplitude  of  the 
errors  only  changed  the  deviation  from  the  nominal  730  km  apogee, 
but  did  not  alter  the  shape  of  the  figures. 

Random  Error 

Randomly  generated  errors  in  navigational  information  were 
investigated  to  determine  if  there  wan  a significant  difference 
from  the  results  obtained  with  constant  and  time- varying  sinusoidal 
errors,  tt  was  also  desired  to  find  out  if  the  correlation  time 
effected  the  results  appreciably.  Correlation  time,  it  will  be 
recalled  from  the  previous  chapter,  is  the  spacing  between  the 
randomly  generated  error  parameters.  Accordingly,  the  program  wa3 
modified  to  run  10  randomly  generated  velocity  error  profiles 
(uniformly  distributed  between  +20  m/sec)  at  each  of  10  correlation 
times.  The  results  are  shown  in  Fig.  12.  The  correlation  time  is 
plotted  versus  the  arithmetic  mean  and  the  standard  deviation  of 
the  apogee  deviations  from  730  km  for  the  10  tra.leetori  es. 

Tt  can  be  seen  that  the  standard  deviation  of  random  error  is 
less  than  the  apogee  deviation  due  to  constant  errors,  again 
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indicating  that  plus  and  minus  errors  tend  to  cancel  each  other  out. 

As  miyht  be  expected,  the  random  errors  produce  greater  apogee 
deviations  than  the  sinusoidal  errors. 

The  relationship  of  correlation  time  to  the  standard  deviation 
is  ambiyuoua.  Although  the  standard  deviation  seems  to  increase 
with  correlation  time  (as  did  the  apoyee  deviations  with  increasing 
period  length  in  the  .sinusoidal  case),  there  is  a minimum  at  25 
seconds.  At  any  rate,  there  does  . ot  seem  to  be  a strong  relation- 
ship between  correlation  time  and  apogee  error. 

Error  Counllnr 

So  far,  in  the  course  of  this  study,  errors  in  velocity,  flight 
path  anyle,  and  acceleration  have  been  input  singly  to  the  simulation. 
It  was  desired  to  determine  if  it  would  be  possible  to  input  combin- 
ations of  the  error  parameters  and  predict  the  apoyee  deviation. 

Usiny  the  principle  of  superposition,  it  would  he  possible  to 
determine  the  deviation  for  each  individual  error,  then  sum  them  to 
find  the  total  apoyee  deviation.  Unfortunately,  this  did  not  prove 
to  be  the  case.  The  results  of  several  tests  indicate  that  the  errors 
are  strongly  couoled  and  superposition  does  not  hold.  This  can  be 
seen  from  the  sample  trajectory  of  Appendix  C.  Constant  errors  of 
-10  m/sec,  .002  radians,  and  .1  m/sec2  were  innut  to  the  simulation. 
From  Figures  6,  7,  and  8,  it  can  be  determined  that  the  apoyee 
deviations  for  each  case  are,  respectively,  +40  km,  -15  km,  and  +12  km, 
totalin"  +87  km  apoyee  error.  However,  the  actual  apoyee  altitude  is 
727  km,  or  -8  km  error.  Mo  further  investiyation  of  the  error 
coupliny  was  attempted,  althourh  the  interrelationship  of  error 
parameters  can  be  seen  from  the  partial  derivatives  of  eq,  (20). 


VT.  Conclusions  and  Recommendations 


Conclusions 

This  study  has  shown  to  what  accuracy  the  AMOOS  vehicle 
navigation  system  must  provide  velocity,  flight  path  apple,  and 
acceleration  information  to  successfully  implement  the  puidance 
scheme.  Thi3  information  can  be  used  to  select  and  desipn  the 
navigation  sensors  for  the  vehicle. 

Velocity  information  must  be  supplied  to  an  aocuraey  of  ±22  m/sec. 
Constant  errors  in  velocity  relate  to  apopee  deviations  aecordinp  to 
the  sensitivity  coefficient  of  -4,25  km/m/sec.  Flipht  path  an pie 
information  must  be  supplied  within  the  accuracy  ranpe  of  -.26  dep 
to  +.20  dep.  Errors  in  measurinp  acceleration  alonp  the  flipht  path 
(which  is  used  to  determine  atmospheric  density  and  then  a density 
altitude)  must  bo  within  + .6  m/seo2.  Errors  in  apopee  altitude 
are  related  to  errors  in  sensinp  acceleration  by  the  approximate 
sensitivity  coefficient  of  1.20  km/m/nee2, 

These  accuracy  tolerances  do  not  seem  to  be  severe  in  view  of 
the  present  capability  of  space  navipation  systems.  The  exoatmos- 
pheric  navipation  problem  imposes  more  strinpent,  accuracy  reouire- 
ments  on  navipational  hardware,  so  there  should  not  be  a major 
problem  in  ohtaininr  the  accuracy  required  for  the  atmospheric  flipht. 
The  puidance  scheme  has  demonstrated  its  capability  to  handle 
time-varyi nr  errors,  both  sinusoid a1  and  random.  There  is  no 
frequency  instability  in  the  puidance  equations,  and  the  anoree 
deviations  for  error  freauencies  above  .1  rad/sec  are  minimal. 
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Recommendations 

The  next  step  in  desifrninr  the  navigation  and  Guidance  system 
for  the  AMOOS  vehicle  would  be  to  select  navigational  sensors 
which  can  provide  the  required  accuracy,  both  for  the  exoatmos- 
pheric  and  the  atmospherio  phases  of  the  mission.  Information  from 
these  sensors  would  have  to  be  processed  by  the  navigation  computer 
before  beinp:  input  to  the  pruidance  system. 

Rxistinm  hardware  would  have  to  be  evaluated  to  determine  i.f 
it  would  be  suitable  for  use  on  the  vehicle,  or  if  new  hardware 
would  have  to  be  designed. 
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APPENDIX  A 


The  AMOOS  Conce pt 


The  Aeromaneuvering  Orbit- to-Orbit  Shuttle  (AMOOS)  is  designed 
to  be  carried  to  and  from  low  earth  orbit  inside  the  Shuttle  Orbiter 
payload  bay.  Basically,  AMOOS  differs  from  the  Space  Tup  in  that  it 
returns  to  low  earth  orbit  by  making  a braking  pass  through  the 
atmosphere.  The  added  vehicle  weight  due  to  increased  structural 
strength  and  the  addition  of  a thermal  protection  system  is  more 
than  offset  by  the  weight  of  fuel  required  for  the  same  propulsive 
braking.  Thus,  AMOOS  can  carry  2 to  3 times  the  payload  to  GKO  for 
the  same  weight  vehicle. 

The  AMOOS  concept  has  been  proposed  for  both  manned  and  un- 
manned missions.  The  mission  of  primary  interest  is  to  geosynchron- 
ous altitude,  and  that  is  the  mission  profile  examined  in  this  study. 
However,  the  concent  is  not  limited  to  high  earth  orbits.  An  AMOOS 
vehicle  could  be  used  upon  returning  to  earth  from  lunar  missions,  or 
it  could  be  used  to  provide  aeromaneuvering  in  any  planetary 
atmosphere, 

In  order  to  provide  for  the  maximum  efficient  use  of  the  Shuttle, 
AMOOS  is  sized  tn  the  maximum  Shuttle  payload  weight.  The  dimensions 
of  the  payload  bav  require  that  AMOOS  be  shared  like  a cylinder  with 
one  end  modified  to  an  aerodynamic  nose  can.  The  general  configuration 
is  shown  in  Fig,  13. 

The  nose  cap  is  hinged  to  permit  the  rocket  engine  to  fire  forward, 
This  aft  cargo  bay  arrangement  permits  the  vehicle  to  be  adapted  to 
modular  configuration.  A hinged  flan  at  the  rear  provides  the 

SQ 


required  aerodynamic  trim  characteristics.  An  ablative  thermal 
protection  system  is  used  on  the  vehicle,  and  it  must  be  replaced 
after  each  mission. 

The  possibility  exists  of  plaeinm  heavy  payloads  into  SEO 
by  staffing  2 AMOOS  vehicles  (usinm  2 Shuttle  flights).  This  would 
permit  the  placing  of  a 15(000  km  payload  into  OBO.  A single  stame 
AMOOS  can  carry  a round  trip  payload  of  2fi60  km,  or  deliver  5200  km. 
This  compares  to  1100  km  and  3600  km,  respectively,  for  the  all- 
propulsive  Space  Turn.  Thus  it  can  be  seen  that  AMOOS  enjoys  a 
definite  performance  advantame  over  the  Space  Turn. 

Estimated  oosts  (Ref  Is 54)  for  the  riesimn,  developement , 
testinm,  evaluation,  and  production  of  AMOOS  are  as  follows  (in  1970 
United  States  dollars); 


nmvMs 


5H.  million 


1st  unit  production  23. 7 million 


refurbishment 
(per  flimht) 


1.1  million 


These  estimates  do  not  include  enmine  costs,  estimated  at  $130 
million  for  DDTAR,  and  $.7  million  unit  cost;  facility  and  operations 
costs;  and  the  prime  contractor's  fee. 

Bimure  14  shows  the  orientation  of  the  vehicle  with  resnect 
to  anmle  of  attack  (d),  bank  anmle  (R),  and  flimht  path  anmle  (v). 

Note  that  hank  anmle,  as  defined,  is  not  the  anmle  produced  by  a 
simple  roll  about  the  vehicle  centerline.  The  autopilot  must  orient 
the  vehicle  to  ohanme  the  direction  of  the  lift  vector  by  anmle  R 
while  keepinm  the  lift  vector,  vehicle  centerline,  and  the  relative 
wind  vector  in  the  same  plane.  Phis  allows  the  vehicle  to  keep  the  same 
side  (covered  with  ablative  material)  flyinm  into  the  relative  wind. 


Definition  of  Angle  of  Attack.  Bank  Angle,  and  Flight 
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Program  Listing 

The  computer  listing  which  follows  was  used  to  simulate  the 
vehicle  trajectories.  The  program  uses  the  following  input 
information  s 

1.  ERRVET,  - velocity  error  » in  this  example,  -10  m/sec  constant 

error 

2.  ERROAM  - flight  path  angle  error  " in  this  example,  +.002  radians 

constant  error 

3.  ERRACC  - acceleration  error  - in  this  example,  +.1  m/sec^ 

constant  error 

4.  TFOR  - format  parameter  - 0 for  short  output,  1 for  full  output 
1,  TTRA.I  - trajectory  designator  - in  this  example,  #1 

6.  VAP  - vacuum  perigee  - 71.1133  km 

7.  DX  - guidance  cycle  time  - 4 sec 

8.  DT  - simulation  time  step  • 1 sec 

9.  C\(t)  - feedback  gain  vector  ( see  Fig.  4) 

The  output  from  this  program  is  given  in  Appendix  C. 
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The  sample  trajectory  contained  in  this  section  was  produced 
b”  the  program  listed  in  Appendix  R. 

The  output  consists  of  the  following: 


1.  Echo  print  of  input  data,  plus  calculated  entry  conditions, 


2.  Nominal  trajectory  whioh  the  cuidance  system  will  try  to  control 
the  vehicle  alone.  Total  time  of  flirht  is  420  seconds,  and 
perigee  (when  flight  path  ancle  chances  from  - to  +)  is  70.113  km. 

3.  Desired  aoocee  altitude  and  time  required  to  descend  to  95.4  km. 

4.  Echo  print  of  feedback  rrain  vector. 

5.  Transfer  orbit  elements  and  inertial  position  and  velocity 
at  becinninc  of  atmospheric  phase  of  flicht. 

6.  Echo  print  of  error  inputs. 

7.  Actual  vehicle  trajeotory,  civinc  actual  altitude,  velocity,  and 
flight  path  ancle;  calculated  densitv;  commanded  bank  ancle; 
and  inclination  of  the  osculatinc  orbit, 

8.  Density  profile  whioh  was  used  in  this  traieotory,  and  the 
actual,  apocee  altitude  the  vehicle  reached. 
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Robert  -1.  Chambers,  .ir.,  was  born  on  Mav  30th,  1948  in  Oranre, 

Mew  Jersey.  He  graduated  from  Mountain  Hirh  School,  West  Oranre,  New 
Jersey  in  June,  19 66  and  entered  the  United  States  Air  Force  Academy, 
Colorado.  He  graduated  in  June,  1970,  receiving  a commission  as  a 
Second  Lieutenant  and  a Haehelor  of  Science  Decree,  with  majors  in 
Astronautics  and  in  Fnsineerinr  Sciences, 

After  attending  pilot  training  at  Crair  AFH,  Alabama  he  was 
assigned  to  Rhein-Main  Air  Rase,  Federal  Republic  of  Oermany. 

While  stationed  here  he  flew  as  a squadron  pilot  in  the  C-130  aircraft, 
and  later  the  C-13l/l'-29  aircraft.  Tn  1973  he  was  selected  to 
attend  the  Air  Force  Institute  of  Technology  and  entered  in  Aurust 
of  that  year. 
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20.  ABSTRACT  (CrniNnud  on  ravaraa  *lda  II  nacaeamty  and  Idanllfy  by  block  numbar) 

" '"A  preliminary  investigation  was  conducted  to  determine  the  navigational 

accuracy  required  by  the  Aeromaneuverinr  Orbit- to-Orbit  Shuttle  (AMOOS) 
during  the  atmospheric  nhase  of  flimht.  The  guidance  scheme,  which  is  the 
same  as  the  one  developed  by  lockheed  in  the  original  A MOOS  study,  uses  the 
parameters  of  velocity,  flight  path  anrla,  and  density  altitude  to  correct 
to  a nominal  tra.iectory.  Density  altitude  is  obtained  from  atmospheric 
density,  which  is  calculated  from  vehicle  acceleration.  Errors  in  velocity, 
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flight  path  angle,  and  acceleration  were  introduced  into  a three-degree 
of  freedom  computer  simulation  of  the  vehicle  trajectory  using  bank 
angle  oommands  generated  by  the  guidance  equations.  Deviations  from 
standard  atmospheric  density  were  taken  into  account.  ' . 4 

The  maximum  errors  allowable  that  still  permitted  tlie  vehicle 
to  attain  its  phasing  orbit  apogee  altitude  of  720  km  WOO  km)  were 
determined.  Three  types  of  error  were  investigated;  constant, 
sinusoidal,  and  random.  For  time- varying  errors  the  frequency  dependence 
was  examined,  as  was  cross-coupling  of  errors. 

It  was  concluded  that  the  guidance  scheme  can  tolerate  fairly  large 
errors  and  still  guide  to  an  acceptable  apogee  altitude.  The  navigational 
accuracy  required  for  the  atmospheric  phase  of  flight  is  within  the 
capability  of  present  day  astrionics. 
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